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Abstract. - We study the nature of the double-peak specific heat structure in kagome clusters. 
That containing 12 spins is considered thoroughly by numerical diagonalization. Simple models 
are proposed revealing the low-T peak nature at T% < A (A is the spin gap) in this case and 
in those of larger clusters studied so far. We show that the rapid increase in density of states 
just above the spin gap gives rise to this peak. These models establish the reason for the weak 
magnetic field sensitivity of the low-T peak. Our discussion could be appropriate for other 
frustrated antiferromagnetic systems too. 



Introduction. — Much attention in last decade was given to such frustrated magnets as 
kagome and pyrochlore. This interest was motivated by their unusual low-T magnetic proper- 
ties observed experimentally (see [21 113) and references therein). For instance, the deviation of 
their uniform susceptibility from Curie- Weis law occurs at much lower temperatures than for 
non- frustrated ones. Then the specific heat C measurements in spin-| kagome material SrCr- 
GaO revealed a peak at T w 5 K which was practically independent of magnetic field up to 12 
T and CocT 2 at T < 5 K £Q. For qualitative understanding of these systems low-T physics 
numerous finite clusters diagonalization studies were carried out 0E1E1EE1EHZH11E1E1 • 

In particular, investigations of kagome Heisenberg spin-i clusters with the number of sites 
N < 36 and with periodic boundary conditions 0J|Sl|ljJ[3IBl revealed a gap A separating 
the ground state from the upper triplet levels and a band of nonmagnetic singlet excitations 
inside the spin gap. It was obtained that the number of states in the singlet band increases 
exponentially with the number of sites. The similar picture with many singlets inside the spin 
gap exists in pyrochlore clusters too 0. 

Specific heat calculations of kagome clusters with N = 12, 18, 24 and 36 revealed two 
peaks with the low-T one at T; < A 001113111] It was obtained for N = 18. 36 in [TT] 
that the low-T peak is weakly dependent on magnetic field and C cx T 2 in a small interval 
below it. The point of view is accepted now that the wealth of low- lying singlet excitations is 
responsible for the low-T peak |10llll| and for its field independence [Tllllj . At the same time 
it was pointed out in [TT1I9] that upper triplet levels contribute to it as well. 

The double-peak specific heat structure is a general feature of many frustrated antiferro- 
magnetic systems. It was obtained numerically in Heisenberg spin-i pyrochlore slab |12j. in 
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A chain and in clusters of triangular lattice in the model with multi-spins exchange taken 
into account |15U14) . The last approach was proposed to describe the second layer of 3 He 
absorbed on graphite [T4] . 

In the present paper we study the nature of the double-peak specific heat structure in 
Heisenberg spin-^ kagome antiferromagnetic clusters. Our conclusions on the origin of the 
low-T peak contradict to those having been proposed in the previous works [10111 1| . We start 
with numerical diagonalization of the kagome Heisenberg spin-^ cluster plotted in the inset 
of Fig. U (star). As is shown in our recent papers the st ar nas doubly degenerate singlet 
ground state separated from the lower triplet level by the gap A. The lower singlet band 
obtained in the previous finite clusters studies 0JE1E10IB1 appears in this approach from the 
stars degenerate ground states as a result of inter-star interaction. 

So there are no singlet states inside the spin gap in the star as it was in clusters with 
periodic boundary conditions studied so far. Nevertheless we have obtained double-peak 
structure of C with the low-T peak at T < A which possesses weak field sensitivity similar 
to results of the previous studies 8,9, 10, 11 . Simple models revealing the low-T peak nature 
is proposed for the kagome star and other larger antiferromagnetic clusters. It is shown that 
the rapid increase in density of states just above the gap is responsible for this peak. These 
models establish the reason for the weak magnetic field sensitivity of the low-T peak. 

The nature of the double-peak structure in other frustrated antiferromagnetic systems can 
be the same as it is for kagome clusters considered here. We hope that the present paper will 
stimulate the corresponding studies. 

Analysis of the star. — We begin with consideration of the Heisenberg spin-^ antiferromag- 
netic cluster shown in the inset of Fig. ^ which Hamiltonian has the form 

h x s < s - "Y. s -- w 

where denote nearest neighbors and the value of exchange J = 1 so as temperature T 
and magnetic field H to be measured in coupling constant unite. We start with H = in 
eq. Q). As the Hamiltonian commutes with all projections of the total spin operator, all star 
levels are classified by the values S, irreducible representations (IRs) of its symmetry group 
Cq v and are degenerated with S z . In the basis of IRs the matrix of the Hamiltonian has 
a block structure. Each block has been diagonalized numerically. Some low-lying levels are 
presented in Table [I] 

As is seen the star has doubly degenerate singlet ground state separated from the lower 
triplet level by the spin gap A w 0.26 ^H]. Energies of the ground state and the upper level 
are —4.5 and 4.5, respectively. A distinguishing characteristic of the spectrum in the range 
A — 4.5 < E < 4.5 is that levels are very close to each other: distances between them are of 
the order of 0.1 -j- 0.01 < A. A part of the spectrum shown in Table [I] reflects this feature. 
We show below that this peculiarity plays the crucial role for low-T star properties. 

Specific heat C calculated with this spectrum is shown in Fig. ^ R nas a double-peak 
structure with low-T peak at T; w 0.085 and high-T one at Th « 0.6. The most intriguing 
feature of C is the existence of the low-T peak at T; < A. The matter is that according to the 
common point of view specific heat of a fully gaped system should decrease exponentially with 
T in the range T < A. Meanwhile in our case the exponential decay occurs at T < T; only. 
Moreover in a usual situation a peak in C results from a peak in the DS. As is demonstrated 
below, it is not the case for the star. We show now that the extremely high DS just above the 
gap is in the origin of this surprising behavior. 

A graphic of integrated density of states (IDS) N(E) which is the number of states with 
energies lower than E is presented in Fig. [21 It is well fitted by the function Af(E) = 
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Table I — Low-lying levels 
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2200tanh(0.39£T) + 2060. As distances between levels at A - 4.5 < E < 4.5 are of the 
order of 0.1 -j- 0.01 one can use this function as IDS above the gap for the specific heat mod- 
eling in the range T ^> 0.01. Then the DS given by p(E) = dN(E)/dE can be approximated 
on the interval A — 4.5 < E < 4.5 as follows: 




Fig. 1 - Specific heat C(T) per spin (solid line) and entropy S(T) divided by <S(oo) = AT In 2 (dashed 
dotted line) of the antiferromagnetic cluster with N = 12 shown below the curves (star). The specific 
heat calculated with eq. iJ^J is shown by dashed line. 
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Fig. 2 - Integrated density of states of the star M(E) which is the number of states with energies 
lower than E. It is fitted by 2200 tanh(0.39i5) + 2060. The inset shows the low-energy sector where 
M{E) is fitted by 180(£ + 4.24). 



where w = 858 and g = 0.39. In this model C has the form 

c _dE d ( dEe~ E / T Ep(E — 4.5) \ 
dT dT\ d+ Jl dEe~ E l T p(E - 4.5) )' 

where p(E) is given by eq. @ and d is degeneracy of the ground state which is equal to 2 
for the star. In eq. and below we shift the energy scale for convenience so as the ground 
state energy to be equal to 0. Specific heat calculated from eqs. @ and © by numerical 
integration is also presented in Fig. ^ It nas a double-peak structure and coincides well 
with that obtained with the real spectrum. It reproduces the high-T peak which evidently 
corresponds to the maximum of the DS at E — and has the low-T peak at T; = 0.1 which 
is 15% higher than the real one. In any way the result reproduces the characteristic features 
of the specific heat and encouraged us to use this model in further C analysis. 

Low-T peak. — Let us consider in detail the specific heat at T < A. We derive now an 
analytic expression for it in this case in the assumption ^ 3> g which is held for the star 
(I «4 and g ss 0.4). It is easy to show using eq. (J2J that at these conditions one can replace 
p(E — 4.5) in eq. © by a constant W w p(A — 4.5) « 117. Unfortunately in this low-energy 
region approximation eq. J2J) is not very good. As is shown in Fig. |3 an accurate value of W 
is 180 because the real Af(E) is well fitted by 180(E + 4.24) in the range A - 4.5 < E < -2.7. 
As a result we have for the specific heat at T < A from eq. J3J) : 



C = 

where x = A/T and Z = d/(WA). 
dC/dx — which has the form 



1 + Zxe x (x 2 + 2x + 2) 



(1 + Zxe x ) 2 



(4) 



Position of the peak 1} is determined by the condition 



Ze x 



(5) 



As follows from eq. © the peak appears if Z < Z c « 2. For the star Z ss 0.04 and numerical 
solution of eq. 10) gives Ti = 0.093 which is very close to the real value of 0.085. Analysis of 
eqs. and JSJ shows that the peak becomes smaller and T\ tends to zero as A decreases. At 
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the same time T becomes larger as W decreases and even exceeds the value of A when W is 
small enough. In any way the peak disappears as soon as Z amounts to Z c . So the value of 
A and that of the jump in the DS above the gap are responsible for the peculiar low-T star 
properties. 

Field dependence. — We turn now to the discussion of the low-T peak dependence on the 
external magnetic field H . In an isotropic unfrustrated antiferromagnet a singularity in C /T 
implies transition to long-range order. It is suppressed to zero temperature by H of the value 
of the peak temperature. In this respect the field behavior of the star at small T is also 
unusual. The results of C calculations using the spectrum of the Hamiltonian eq. obtained 
numerically are presented in the left graph of Fig. [3J It is seen that at H = Ti « 0.085 the 
peak height is reduced 11% only and T is diminished slightly. The field of the value of the 
spin gap is needed to wash out the peak. 

The low-T peak evolution can be described by the model discussed above. It allows also 
to understand the nature of this weak H sensitivity. The field splits levels with S ^ 0. It is 
clear from Table [i] that the low-energy sector above the gap is represented mostly by triplets. 
So the DS can be modeled at H < A in the following way: it is at < E < A — H, the 
DS is W/3 at A - H < E < A, it is 2W/3 a,tA<E<A + H and the DS remains W 
at E > A + H. The results of the peak evolution with H in this approach are presented in 
the right graph of Fig. |3| It is seen that this model reproduces even quantitatively the main 
tendencies in low-T peak evolution with H . The field smooths the area where the drop in DS 
takes place and in doing so it causes the peak reduction. For the changes in the peak to be 
significant H « A is needed. So the reason for the weak field sensitivity of the low-T peak in 
the star is that T; w A/3. 

Larger kagome clusters. — As was mentioned above the double-peak specific heat structure 
in kagome clusters with even N was obtained in many previous numerical works [51819111011 lj . 
Height and position with Th~2/2, of the high-T peak got by high-temperature expansion ^H] 
and by finite cluster diagonalization studies |5ll%ll§ HT(lllll| coincide with each other and with 
our results. The low-T peak with T; < A was obtained in [HHHIEII^J but its position and 
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height were depended on N, cluster form and calculation technique. We show now that the 
reason for the low-T peak appearance in those cases is similar to that for the star discussed 
above: rapid increase in the DS. It will be argued that the low field sensitivity of the low-T 
peak obtained in has the same origin as it is for the star. 

We focus here on the largest cluster with N = 36 studied before which low-energy spectrum 
was obtained in [7||^j and thermodynamical properties were discussed in According to p[] 
the spectrum above A is nearly continuous as it is for the star. But in contrast the spin gap 
there is filled with singlets separated from the non-degenerate singlet ground state by a very 
small gap A s <C A. According to results of [7] distribution of these singlets is quite even with 
DS of the order of 1000. The DS above the gap is much larger than that inside it. 

Qualitative understanding of the low-T peak nature in this case can be obtained by the 
model similar to that proposed above for the star. Let us assume that the DS is W at E < A', 
where A' > A, and there is a jump in the DS at E = A' so as it is equal to W ^> W above 
A'. Specific heat calculation based on the formula similar to eq. © gives on the interval of 
interest A s « T < A' the following equation: 

c _(l + re x f + rx 2 e x + 2rZxe 2x + Zxe x (x 2 + 2x + 2) 

_ (1 + re x + Zxe x ) 2 ' ^ ' 

where now x — A'/T, Z = l/(A'W) and r — W'/W. We have made an effort to reproduce 
the low-T peak in C/N obtained in JI] with the height of approximately 0.152 and T ~ 
0.05 < A w 0.074. eq. © at A' = 0.23, W = 43000 and W = 1000 gives a peak with 
these parameters but the approach is too rough to get T 2 behavior of C at T < T. As is 
demonstrated below there are no more sets of parameters in this model describing the peak 
discussed in In this case ZA'/Tj -C r and terms in eq. © containing Z do not affect the 
low-T peak. As a result T\ is given by the equation 

re x = — !— . (7) 
x-1 y ' 

So the low-T peak properties are determined by A' and the ratio r. Thus they are a result 
of cooperative effect of states above and below A'. This finding is in accordance with that 

of m. 

It is seen from eqs. (@ and (7J that value of A' does not affect the peak height but it 
determines T/. Evidently T decreases as A' tends to zero. The value of r determines both the 
peak position and its height. The peak becomes smaller and T increases as r becomes larger. 
So it is clear now that if, according to [7J, DS below A' is chosen equal to 1000 and W <C W 
the values of A' and W are determined unambiguously by the low-T peak obtained in |11|. 

Using this model it is easy to describe qualitatively the weak sensitivity of C to the 
magnetic field which was discussed in JT]. It was demonstrated there that at H = Ti = 0.05 
the peak is decreased by 10% only and T is also decreased slightly. The situation here is 
similar to that for the star considered above. As magnetic field splits levels with S ^ the 
DS in the field H < A' will be the following: the DS at < E < A' - H remains W, it is 
W + W/3 at A' - H < E < A', the DS is 2W/3 at A' < E < A' + H and the DS remains 
W above A' + H . We also assume here that the spectrum above A' at H — is formed 
mostly by triplets. As a result of C calculations in our model it was obtained that the field 
H = Ti = 0.05 reduces the low-T peak by 8% and diminishes T slightly. These findings are 
in even quantitative egreement with those of The reduction of the peak at H = A'/2 
calculated in our model is about 28%. So we see that field H ~ A' ^> T is needed to change 
the low-T peak significantly. 
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Entropy analysis. — Entropy of the star S is also presented in Fig. ^ Its temperature 
dependence obeys the most characteristic feature of those of the larger clusters: in all cases 
there is 50% of the total entropy in the low-T peak at T < 0.2 [HUTU] . It is seen from the 
above consideration that this pecularity in kagome clusters stems from the rapid increas in 
DS and states above A' give the main contribution to S at T < 0.2. Such an entropy feature 
was obtained in other models of frustrated antiferromagnets mentioned above and having 
double-peak C structure [T31IT1) . 

In closing we point out that the low-T peak in C can remain in spin-^ kagome anti- 
ferromagnets. The peak at small T with weak field sensitivity obtained experimentally in 
SrCrGaO could be of the same nature. 

Conclusion. — In this paper we study the nature of the double-peak specific heat structure 
in kagome clusters. That containing 12 spins and shown in the inset of Fig. ^ (star) is 
considered thoroughly by numerical diagonalization. Simple models are proposed revealing 
the low-T peak nature at Tj < A (A is the spin gap) in this case and in those of larger clusters 
studied so far numerically (HJEIEIIE]- We show that the rapid increase in density of states 
just above the gap gives rise to this peak. These models establish the reason for the weak 
magnetic field sensitivity of the low-T peak. Our consideration could be appropriate for other 
frustrated antiferromagnetic systems too and, as we hope, will stimulate further studies. 
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